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A procedure to predict the occurrence of periodic clnsters in a system of globally coupled maps displaying a constant mean 
field is presented. The method employs the analogy between a system of globally conpled maps and a single map driven by 
a constant force. By obtaining the asymptotic orbits of the driven map, an associated conpling function can be constructed. 
This function allows to establish a direct connection between both systems. Some applications are shown. 
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The phenomenon of emergence of nontrivial collective 
behaviors in dynamical systems of interacting chaotic el¬ 
ements has been a focus of recent attention [Kaneko, 
1990; Chate & Manneville, 1992(a), 1992(b); Chate et 
aL, 1996; Chawanyafe Morita, 1998]. An important class 
of systems which can exhibit ordered collective effects 
are globally coupled nonlinear oscillators. Such systems 
arise naturally in the description of Josephson junctions 
arrays, charge density waves, multimode lasers, neural 
dynamics, ecological and evolution models [Wiesenfeld 
& Hadley, 1989; Wiesenfeld et al., 1990; Strogatz et al, 
(1989); Nakagawa & Kuramoto, 1994; Kaneko, 1994]. 
Many of these systems can be modeled with globally cou¬ 
pled maps (GCM) [Kaneko, 1990]. In particular, GGM’s 
can show a variety of collective behaviors such as: a) 
formation of clusters, i.e., differentiated subsets of syn¬ 
chronized elements within the network [Kaneko, 1995]; 
b) non-statistical fluctuations of the mean held of the en¬ 
semble [Kaneko, 1995; Perez et al., 1993]; and c) global 
quasiperiodic motion ]Kaneko, 1991; Pikovsky & Kurths, 
1994]. 

Usually, numerical simulations have been performed on 
the entire globally coupled system to search for the oc¬ 
currence of various types of ordered collective behaviors. 
However, this direct procedure gives little information 
about the mechanism for the emergence of collective be¬ 
haviors. In a previous work [Parravano & Gosenza, 1998] 
we have proposed an alternative method that allows to 
predict the emergence of given types of periodic collective 
behaviors in GGM systems exhibiting periodic clusters. 
In this article, we investigate the conditions for the emer¬ 
gence of a class of ordered collective behavior in a GGM 
characterized by the occurrence of K identical clusters, 
all of them runing periodically among K states but out 
of phase with respect to each other in such a way that the 
mean field of the system remains constant during its evo¬ 
lution. This kind of collective behavior has been observed 
for globally coupled logistic maps as three identical out 
of phase clusters by Shinbrot [1994], who noticed that 
the addition of a constant term to a single local map has 
a similar dynamics as the GCM system. Here we use the 


similarity between a steadly driven map and a GCM to 
predict the conditions for the occurrence of general clus¬ 
tered collective behavior in GCM’s displaying a constant 
mean field. 

We consider the globally coupled map system 

= (1 - e)f{x\) + e{xt ), (1) 

with mean field coupling 
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where xl (i = 1,..., N) gives the state of the fth element 
at discrete time t; N is the size of the system; e measures 
the strength of the coupling; and f{x) prescribes the local 
dynamics. Note that the mean held coupling affects all 
the elements in (Eq. ®) in exactly the same way at 
all times. This property of the coupling term allows to 
establish an analogy between the GCM system [Eq. ®] 
and an associated driven map 

St+i = (1 “ £)/('St) + eLt , (3) 

where St is the state of the driven map at discrete time 
t, f{st) is the same local dynamics as in Eq. ®, and Lt 
is the driving term. In general, L* can be any function 
of the discrete time t. Then, any element xl in the GCM 
system behaves equivalently to a single driven map if 
Lt = (xt), Vt, and their initial conditions are identical, 
i.e., xl = So. As it will be shown below, the analogy can 
also be drawn if Lt = (xt) only for asymptotic times. The 
basic idea is that, by studying the asymptotic dynamics 
of the associated driven map, one may infer if specihc 
types of ordered collective behaviors can emerge in the 
GCM. 

The simplest situation for which the associated driven 
map [Eq. (^] can be used to predict the emergence of or¬ 
dered collective behavior in the GCM system [Eq. ®] oc¬ 
curs when the mean field remains constant, i.e., (xt) = C. 
This behavior takes place in the GCM when K identical 
clusters, each having N/K elements and period K, are 
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evolving out of phase in order to yield a constant (xt). 
In this case, the behavior of any of such clusters in the 
GCM, can be emulated by an associated driven map sub¬ 
jected to a constant force L* = C. 

The steadly driven map will have a unique asymptotic 
orbit whenever the local map / satisfies Singer’s theo¬ 
rem. This is the case for the logistic map f{x) = 1 — ax^, 
a G (0,2], which will be used as an example to explain 
our procedure. The unique asymptotic orbit of the driven 
map can be periodic on some regions of its space of pa¬ 
rameters (i.e., in (a,e, C) for the chosen /). Let us con¬ 
sider an orbit of period K, and let the sequence of the 
points on this periodic orbit be {cti, tT 2 , ■ ■ •, o-r}- By us¬ 
ing these K points on the periodic orbit of the driven 
map, the following associated coupling function can be 
defined 
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0(a,e,C') = — . (4) 

i=i 

The associated function 0 depends on the values 
{o’!, ( 72 , ■ ■ •, ixk}, which themselves depend on the value 
of the constant drive C, on the coupling parameter e, 
and on the local map f{x). The function 0 gives the 
value that the mean field in the GCM system [Eq. ( 0 )] 
would have at a given time if at this instant its elements 
are segregated in K identical clusters distributed in the 
states {(Ti,..., (Tif}. Note that the relation Q = (xt) will 
persist for subsequent times if the clusters in the GCM 
sequentially follow the values {(Ti, ..., ctk} while remain¬ 
ing shifted one time step with respect to each other. The 
relation (xt) = 0 is guaranteed for all times when the 
driven map satisfies the condition 0 = C. If C* is the 
solution of 0(a, e, C) = C, then the evolution of a set of 
K independent maps driven by the same constant force 
C* describes a collective behavior in the corresponding 
GCM consisting of K identical, period K clusters which 
maintains {xt) = The driven maps, as well as the 
clusters in the GCM, will cyclically commute over the 
states {(Ti(C*),..., (7k{C*)}. 

It should be noticed that for a given GCM the associ¬ 
ated coupling function 0 can be constructed only when 
the associated driven map displays periodic asymptotic 
responses for some values of Lt = C. As an example, 
Fig. 1 shows the bifurcation diagram of St as function of 
e for the local map f{x) = 1 — ax"^ with a = 2, driven 
by a constant term C = 0.07. Figure 1 also shows the 
time average s of the iterates St as a function of e. The 
statistical quantity s gives the value of 0 on the periodic 
regions of the bifurcation diagram. The intersection of 
the curve s with the line C = 0.07 in the period three 
window gives the value of e that satisfies the condition 
0(a = 2, e, C = 0.07) = 0.07. Consequently, for this 
value of e (« 0.065) and for a = 2, a GCM system cou¬ 
pled through its mean field can exhibit a collective be¬ 
havior consisting of three identical out of phase clusters 


that yield a constant {xt) = 0.07. 

Note that this prediction is made without direct nu¬ 
merical simulation on the GCM system. However, the 
solution of the equation 0 = C does not tell which initial 
conditions on the GCM will conduce towards this partic¬ 
ular collective behavior. We have verified that the GCM 
system coupled through the mean field reaches the pre¬ 
dicted behavior starting from initial conditions {x),} con¬ 
sisting of three groups of A^/3 elements, each group dis¬ 
tributed around one of the states {ai, ( 72 , ( 73 } with some 
dispersion. For N = 1200, the system reaches the ex¬ 
pected behavior when the dispersion is less than 0.007. 
For larger dispersions or for initially uniform distributed 
elements on the interval [—1,1], the GCM system tends to 
split into three similar but not identical clusters result¬ 
ing in a period three mean field that oscillates around 
the value 0.07 with a small amplitude. The later case 
was observed by Shinbrot [1994] in a GCM system with 
a quadratic local map and mean field coupling. In gen¬ 
eral, oscillations of the mean field arise in the clustered 
collective behavior of a GCM when N is not a multiple 
of AT. 

Figure 2 shows 0(a = 2,e,C) as a function of C for 
several values of the coupling strength e. The continuous 
curves correspond to values of 0 in periodic regions of 
the indicated values of e. In the chaotic regions 0 cannot 
be calculated; however we show s for the three smaller 
values of e for which chaos can occur. The intersections 
of the curves with the diagonal give the solutions to the 
equations 0 (a = 2 , e, C) = C for each e shown. 

Figure 3 shows the solution of the equation 0(a = 
2, e,C) = C on the plane of parameters (e, C) for various 
regions where the associated driven map has periodic or¬ 
bits. The labels on the different portions of the solution 
curve indicate the period K of the orbits. For each value 
of e and C on this curve a collective behavior comprising 
K periodic clusters that result in a constant mean field 
{xt) = C can take place in the corresponding GGM for 
appropriate initial conditions. The inserts in Fig. 3 are 
magnifications of the solution 0 = C corresponding to 
two different periodic windows of the driven map. 

In summary, based on the analogy between a GCM 
system and a single driven map, we have presented a 
method to predict the parameter values for the existence 
of periodic clustered collective behavior in GCM’s dis¬ 
playing a constant mean field. By obtaining the asymp¬ 
totic orbits of the driven map, the associated coupling 
function 0 can be constructed. This function allows 
to establish a direct connection between a GCM and 
its associated driven map. The solutions of equation 
0(a, e,C') = C represent surfaces in the space of param¬ 
eters (a, e, C) of the driven map. For a set of parameter 
values (a, e, C) on these surfaces, the corresponding GCM 
can have a mean field (xt) = C{a, e) resulting from a peri¬ 
odic clustered behavior. The method also gives the states 
{(7i(C),..., (7if (C)} sequentially visited by the clusters. 
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Results were shown for the logistic map with a = 2, how¬ 
ever, the procedure can be applied for other local maps 
and for any global coupling function invariant to permu¬ 
tations of its N arguments. 
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FIGURE CAPTIONS 

Fig. 1 : Bifurcation diagram for Eq. (3) as a function of e, 
with fixed a = 2 and Lt = C = 0.07. The horizontal line 
corresponds to the constant value of C. The continuous 
curve is the time average s, which gives the function 0 
on the periodic regions of the diagram. 

Fig. 2 : Associated coupling function 0 (continuous 
curves) on periodic regions of Eq. (3), as a function of 
C, for several values of e and fixed a = 2. The numbers 
in parenthesis indicate the periodicity K. The dots corre¬ 
spond to values of the time average s on chaotic regions. 

Fig. 3 : The solution of 0(a = 2,e,C) = C for various 
regions where Eq. (3) has periodic orbits. Labels indicate 
the period K of the orbits. The inserts are magnifications 
of the indicated portions of the curve corresponding to 
two different periodic windows of the driven map 
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